Fractal Weyl law for chaotic microcavities: Fresnel's laws imply multifractal scattering 
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I. INTRODUCTION 

Optical microcavities are expected to be key components 
in future photonic applications, such as ultralow threshold 
lasers 01 0], single-photon emitters JH 01, and correlated 
photon-pair emitters [5]. Microdisk cavities with non-circular 
cross-sectional shape have recently attracted considerable in- 
terest in the quantum chaos community since the internal ray 
dynamics is non-integrable Da, 01- In fact, microdisks can be 
considered as open billiards. In a billiard a point-like par- 
ticle moves freely in a two-dimensional domain with elastic 
reflections at the boundary [8]. Depending on the shape of 
the boundary the system shows a variety of dynamical behav- 
iors ranging from integrable to fully chaotic [9]. Light rays in 
a microcavity behave similarly as they are totally reflected at 
the boundary as long as the angle of incidence \ (measured 
from the normal) is larger than the critical angle for total in- 
ternal reflection \c — arcsin (1/n), where n is the index of 
refraction. If, however, \ < \c then the light ray escapes 
refractively according to Snell's and Fresnel's laws. 

An important topic of quantum chaos is the analysis of 
the statistical properties of energy levels in quantum systems 
whose classical dynamics is fully or partially chaotic IToll . 
In the last years, the focus has shifted from closed to open 
systems. Quantum eigenenergies of open systems (resonance 
frequencies in the case of microcavities) are complex valued 
with the imaginary part being related to the lifetime of the 
state. Of particular interest is the fractal Weyl law for open 
chaotic systems d [H [H E El [nil . This is an ex- 
tension of the well-known Weyl's formula for bounded sys- 
tems which states that the number of levels with wave number 
k n < k is asymptotically N(k) ~ k 2 for the particular case of 
a two-dimensional system which scales with the energy such 
as quantum billiards. The fractal Weyl law for an open chaotic 
system (which again scales with the energy) having complex 
wave numbers k n can be written as 

N(k) = {k n : Im(fc„) > -C, Re(fc„) < k} ~ k a . (1) 

The cutoff constant C > removes fast decaying states. It is 



conjectured that the non-integer exponent is 



where D is the fractal dimension of the chaotic repeller of 
the open chaotic system Hill : d = D — 1 is the dimension 
of the repeller in a properly chosen surface of section lfl8ll . 
The chaotic repeller is the set of points in phase space that 
never leave the system both in forward and in backward time 
dynamics, see, e.g., Ref. H. 

An essential prerequisite for a statistical analysis is a suffi- 
cient amount of resonance data. Resonance spectra can either 
be measured experimentally or computed by exact quantum 
or semiclassical methods. For some systems, e.g., atoms in 
external electric and magnetic fields the complex 5-matrix 
poles can be obtained directly by diagonalization of non- 
Hermitian matrices using complex dilated basis sets lfl9ll20ll . 
However, in m any theoretical calculations of, e.g ., open bil- 
liards I21I I22I l23l l24ll . optical microcavities Ipal. or in the 
cycle-expanded Gutzwiller-Voros (^-function |26ll the reso- 
nances cannot be obtained directly by matrix diagonalization 
but rather by a numerically expensive root search 

detM(£ n ) = 0, (3) 

where the matrix 'M.(E) depends on the complex valued en- 
ergy E (or on the complex wave number k for billiards and 
microcavities). The calculation of a real valued spectrum 
a(E) as the superposition of Lorentzian line shapes is much 
less expensive, and it would be a great advantage to be able 
to extract the complex S'-matrix poles from a real spectrum 
ct(E). In principle, this can be achieved by fitting the real val- 
ued spectrum to a sum of Lorentzians. Unfortunately, the es- 
tablished methods for fitting Lorentzians are restricted to iso- 
lated or weakly overlapping resonances, whereas typical open 
systems in physically interesting regions are characterized by 
strongly overlapping resonances. 

In this paper we will apply the harmonic inversion method 
to obtain the resonances of open dielectric microcavities. We 
will show that the conjectured fractal Weyl law for open 
chaotic systems is valid for microcavities only if the concept 
of the chaotic repeller is extended to a multifractal by incor- 
porating Fresnel's laws. The paper is organized as follows. In 
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Sec. |II]we present the numerical technique for the extraction 
of the resonance poles from a real spectrum <r(E). In Sec. Hill 
we introduce the dielectric microstadium. The high-index and 
low-index microstadium are investigated in detail in Secs.lIVI 
andM Conclusion are given in SeclVll 



II. THE HARMONIC INVERSION TECHNIQUE 

Here, we will apply methods for high-resolution signal pro- 
cessing to extract the resonance positions, widths, and ampli- 
tudes from a real spectrum a(E). In a first step the spec- 
trum is Fourier transformed with appropriate frequency filters 
to obtain a band-limited time signal. This time signal is an- 
alyzed, in the second step, by a high-resolution harmonic in- 
version method. The harmonic inversion method for signal 
processing has been introduced by Wall and Neuhauser 112711 
and then refined and improved by Mandelshtam and Taylor 
J28il29tl . It has proved to be a powerful tool for both the high- 
resolution analysis of quantum spectra ll30ll and the semiclas- 
sical quantization of non-integrable systems l3ll [32I l33i l34tl . 
For a review on the use of the harmonic inversion method in 
semiclassical physics see Ref. 13511 . Technically, the harmonic 
inversion problem can be recast as a set of nonlinear equa- 
tions, which can be solved by either linear predictor, Pade 
approximant, or direct signal diagonalization |36]. Although 
details of the procedure have already been published we pro- 
vide here, for the convenience of the reader and to make the 
paper self-contained, a brief description of the harmonic in- 
version method extended to the extraction of resonances from 
experimental or theoretical spectra given as superpositions of 
Lorentzians, i.e., spectra of the form 
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where the {fife, T k , A k } are the frequency positions, widths, 
and amplitudes of the resonances. The frequency spectrum 
(0]l can be interpreted as the Fourier transform of a time signal 
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with complex frequencies oj k = fife — ^k and complex val- 
ues d k simply related to the amplitudes A k . The challenge 
is to determine very broad resonances deep in the complex 
plane and to resolve individual resonance poles in regions with 
strongly overlapping resonances. Basically, the harmonic in- 
version algorithm is split into the following two steps. 

In a first step, a frequency window [luq — Aw /2 , wo + Aw / 2] 
is chosen. A band-limited signal Cbi(t) with a finite num- 
ber of about 50 to 200 frequencies is obtained by the appli- 
cation of a windowed discrete Fourier transform to the fre- 
quency spectrum (@). The band-limited signal is evaluated on 
an equidistant grid with time step r = 2tt/Auj and reads 



with n = 0, 1, ... , A^bi — 1. The limits j\ and of the sum 
in Eq. (O are taken so that the equidistant grid points u>j of 
the digitized spectrum (01 cover the selected frequency filter 
[u>0 — Aw/2, wo + Aw/ 2]. The frequencies in the exponent 
of Eq. © are shifted by — wo to reduce the phase oscillations 
of the band-limited signal. The number iVbi of data points c n 
of the band-limited signal is restricted by the number of grid 
points in the frequency window, i.e., iVbi < % — ji + 1. We 
choose iVbi = 2K, where K is an upper bound of the number 
of resonances in the window (see below). 

In a second step, the (shifted) frequencies wi = u> k — oj of 
the band-limited signal © are obtained by solving the nonlin- 
ear set of equations 
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where z k = exp(— iw fe r) and d k are in general complex vari- 
ational parameters. As the number of frequencies in the band- 
limited signal is relatively small (K ~ 50 — 200), several 
variants of Prony's method |H3l, which otherwise would be 
numerically unstable, can now be applied successfully, e.g., 
linear predictor, Pade approximant, or direct signal diagonal- 
ization J36|]. Here, we briefly elaborate on the Pade approxi- 
mant. 

Let us assume for the moment that the signal points c n are 
known up to infinity, n = 0, 1, . . . oo. Interpreting the c„'s as 
the coefficients of a Maclaurin series in the variable we 
can define the function g(z) = Y^=o c n z ~ n ■ With Eq. © 
and the sum rule for geometric series we obtain 
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The right-hand side of Eq. dHJ is a rational function with poly- 
nomials of degree K in the numerator and denominator. Ev- 
idently, the parameters z k = exp (— iw fc r) are the poles of 
g(z), i.e., the zeros of the polynomial Qk(z)- The parame- 
ters d k are calculated via the residues of the last two terms of 
Eq. ([8j. Application of the residue calculus yields 
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with the prime indicating the derivative d/dz. Of course, the 
assumption that the coefficients c„ are known up to infinity is 
not fulfilled and, therefore, the sum over all c„z~™ in Eq. ^ 
cannot be evaluated in practice. However, the convergence 
of the sum can be accelerated by use of the Pade approxi- 
mant. Indeed, knowledge of 2K signal points cq, . . . , c%k— l 
is sufficient for the calculation of the coefficients of the two 
polynomials 
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The coefficients a k with k = l,...,K are obtained as solu- 
tions of the linear set of equations 

A" 

c n = ^a k c n+k , n = 0, . . .,K- 1 . (11) 

fc=i 

Once the a k are known, the coefficients b k are given by the 
explicit formula 

K-k 

bk = ^ a k+m c m , k = l,...,K. (12) 

m=0 

Note that the Pade approximant does not only accelerate the 
convergence of the sum over all c n z~ n in Eq. <[8J but also 
yields an analytic continuation for z values, where the sum is 
not absolutely convergent. 

The parameters z k = exp(— \oj' k r) and thus the frequencies 

u k = lu + Lu' k = uj + - In z k (13) 

r 

are obtained by searching for the zeros of the polynomial 
Qk(z) in Eq. ( TTOb . This is the only nonlinear step of the al- 
gorithm. The roots of polynomials can be found, in principle, 
by application of Laguerre's method rf38ll . However, it turns 
out that an alternative method, i.e., the diagonalization of the 
Hessenberg matrix 
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with a k the coefficients of the polynomial Qk(z) in Eq. (TTOb 
is a numerically more robust technique for finding the roots of 
high degree (K > 60) polynomials [38]. 

When the harmonic inversion method is applied to a contin- 
uous spectrum of a scattering system, given as a superposition 
of Lorentzian line shapes, the resulting u k are the complex 
poles of the S'-matrix. 



cavity plane. The wave function ip and its normal derivative 
are continuous across the boundary of the cavity. At infinity, 
outgoing wave conditions are imposed. Note that we ignore 
in the following a slight frequency dependence of n. 

There are only a few reports in the literature on the spec- 
tral statistics of microcavities. In a recent work the resonance 
width distribution in the circular disk has been studied l4(ill . 
The integrability of this system allows to easily gather suffi- 
cient resonance data. In Ref. lT4lTl the resonance width distri- 
bution of a circular cavity with strong surface roughness has 
been analyzed. However, an additional approximation was 
necessary in order to get a sufficient amount of resonances. 

In our approach, no such approximation is needed. More- 
over, it is not restricted to a particular geometry. As exam- 
ple, we consider the dielectric microstadium as illustrated in 
Fig-El The cross-sectional area is bounded by two semicircles 
with radius R > and two straight lines of length 2L > 0. 
The corresponding stadium billiard is a paradigm for fully 
chaotic systems l42l 14311 . We choose the case L = R which is 
the most chaotic one 1441. 
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FIG. 1 : Schematic top view of the microstadium cavity. Total inter- 
nal reflection of a ray with angle of incidence \ ls illustrated. The 
arclength coordinate along the circumference of the cavity's bound- 
ary is denoted by s. 

The spatial structure of optical modes in microstadiums has 
been extensively studied in the context of ray-wave correspon- 
dence in open systems ll4^l4^il47ll48ll49ll . Our aim is to study 
the spectral properties of modes in such a kind of cavity. 

IV. HIGH-INDEX MICROSTADIUM 



III. THE MICROSTADIUM 

Below a certain cutoff frequency w cu t ff, microdisk-like 
cavities can be regarded as two-dimensional systems. In this 
case Maxwell's equations reduce to a two-dimensional scalar 
mode equation B9I1 

-S7^ = n 2 (x,y)^ : (15) 
tr 

with effective index of refraction n, frequency ui and the speed 
of light in vacuum c. The mode equation (l5[ is valid for 
both transverse magnetic (TM) and transverse electric (TE) 
polarization. We focus on TM polarization with the electric 
field E(x,y,t) = (0, 0, ip(x, y)e~ lujt ) perpendicular to the 



We consider first a GaAs microstadium with a high refrac- 
tive index n = 3.3, the exp erimental realization of which has 
been reported in Refs. Il50il5lll . 

A. Spectral analysis of resonances 

The modes ([13] ) can be numerically computed with the 
boundary element method (BEM) using a root search algo- 
rithm in the complex frequency plane 12511 . Even though the 
BEM is very efficient, it does not allow for obtaining a suf- 
ficient amount of resonances for a statistical analysis. Our 
improved strategy is therefore based on the BEM and the har- 
monic inversion technique. It consists of two steps. In a first 
step, we compute spectral data on the real frequency axis. For 
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convenience, we choose the total scattering cross section a 
as function of the normalized frequency £1 = u)R/c = kR 
using the BEM. For a definition of the scattering cross sec- 
tion we refer to Ref. l25ll . We are able to cover the inter- 
val fi S [0, 25] where we found 3772 resonances. The upper 
bound = 25 corresponds to R = 3.2 /j,m assuming a typi- 
cal wavelength of 800 nm in GaAs devices Il50l,l5lll . To keep 
the picture clear, Fig. [2] shows the calculated scattering cross 
section in a smaller frequency interval. 




FIG. 2: Calculated total scattering cross section a versus real-valued 
frequency Q, = uiR/c for the high-index microstadium. The solid 
curve shows the result of the BEM with a plane wave incidence at 
30° to the horizontal line segments of the stadium. The dashed curve 
is the reconstructed a based on the resonances extracted by the har- 
monic inversion. 

In a second step we analyze the spectral data on the real 
frequency axis with the help of the harmonic inversion tech- 
nique. Fig. [3] shows the resulting resonances in the complex 
plane with Im(f2) > —0.06. The imaginary part is related to 
the quality factor via Q = — Re(f2)/[2Im(f2)]. The limita- 
tions of the harmonic inversion method are as follows: Very 
weak and broad resonances, i.e., resonances with small ampli- 
tudes Ak and large widths Tk in Eq. (HJl may be overlooked. 
Also, in rare cases one true resonance can be fitted by a sum of 
two Lorentzians with nearly identical resonance positions and 
widths. We have two ways to verify the harmonic inversion 
technique. First, we can check exemplarily resonances with 
the BEM by using root searching in the complex plane. We 
always find good agreement. Second, from the extracted reso- 
nances we can reconstruct the scattering cross section. Fig. [2] 
reveals good agreement with the original scattering cross sec- 
tion. 

The statistics of imaginary parts Im(f2) = —R/(2cr), 
where r is the lifetime of the given mode, is depicted in Fig.|4] 
We see a clear maximum around Im(£l) « —0.014. 

Figure [5] shows the number of modes with cutoff C = 0.06 
according to Eq. (Q]) versus normalized frequency £1 = kR 
in a log-log plot. The data can be extremely well fitted by 
a straight line with slope a w 1.98. Because of the well pro- 
nounced maximum in the statistics of imaginary parts in Fig. [4] 
we expect only a small change of the exponent a when the 
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FIG. 3: Resonances in the plane of complex frequencies SI for the 
microstadium with n = 3.3. We found 3772 resonances in total. 
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FIG. 4: Probability density of imaginary part of resonant frequencies 
Q for the high-index microstadium. 



cutoff C is varied. Indeed, this expectation is confirmed by 
the numerics as we find a £ [1.96, 2.02] for C G [0.03, 0.1]. 
As an exponent larger than 2 for two-dimensional open sys- 
tems is rather unphysical, we believe that the regime [2, 2.02] 
is solely due to the numerical uncertainty. 



B. The chaotic repeller 

The chaotic repeller is the set of points in phase space that 
never leave the system both in forward and in backward time 
dynamics lfl8tl . The stable (unstable) manifold of a chaotic 
repeller is the set of points that converges to the repeller in 
forward (backward) time evolution. The unstable manifold 
therefore describes the route of escape from the chaotic sys- 
tem. For the case of chaotic microcavities it has been demon- 
strated that the unstable manifold ca n p l ay an important role 
for the far- field emission pattern 11481 P2ll53ll54l. 15511 . 

Here, we discuss the chaotic repeller of the microstadium 
restricted to the Poincare surface of section (SOS) lfl8ll . It is a 
plot of the intersection points of a set of ray trajectories with 
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FIG. 5: Number of modes N with cutoff C = 0.06 versus frequency 
SI in a double-logarithmic plot (solid line) for the microstadium with 
n = 3.3. Dashed line with slope a ~ 1.98 is a least-square fit; offset 
for clarity. 

a surface in phase space. For billiards and microcavities it is 
convenient to take the cavity boundary as SOS. The dimen- 
sion of the repeller on this reduced phase space is d = D — 1. 
The reduced phase space is parametrized by the arclength co- 
ordinate along the circumference of the cavity's boundary, s, 
and the tangential momentum, p = sin x, where x is the angle 
of incidence measured from the surface normal; see Fig.Q] 

To compute the repeller we portion the phase space above 
the upper critical line for total internal reflection, p > 
jsinxd = 1/n, into L x L = 1024 x 1024 boxes. For each 
box we would like to determine whether it belongs to the re- 
peller or not. Obviously, this depends on the chosen point 
within the box. To soften this arbitrariness the following aver- 
aging procedure is used. In each box m = 400 trajectories are 
started and propagated both forward and backward in time. 
Given a cutoff length L cuto ff = 20. 5R, we assign an intensity 
J = 1 to rays which never enter the leaky region \p\ < 1/n, 
otherwise 1 = 0. The weight w E [0, 1] associated with each 
box is then 

W = (16) 

i=l 

Figure [6] shows the chaotic repeller determined in this way. 
The magnifications in Fig. [7] demonstrate the fractal nature of 
the repeller. 

In order to compute the fractal dimension of the chaotic 
repeller, we apply the thermodynamic formalism for fractal 
measures; see, e.g., Ref. 15611 . As for the determination of the 
chaotic repeller we divide the region (x, y) = [s/s max , (p — 
l/n)/(l - 1/n)] e [0, 1] x [0, 1] into Lx L= 1024 x 1024 
boxes with size e x e where e = l/L. The weights w^i = 
w(xk, yi) from Eq. ( fTSI l are normalized such that J2 kl Wki — 
1. We write the coarse-grained partition function as 

r(e,g)=5>^, (17) 

kl 
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FIG. 6: Chaotic repeller (dark regions) of the high-index microsta- 
dium in the Poincare surface of section. The horizontal axis is the 
arclength coordinate s/s m ax G [0, 1] and the vertical axis is the tan- 
gential momentum p = sinx £ [1 /n, 1] . The rectangular region is 
magnified in Fig.|7Ja). 




FIG. 7: (a) Magnification of rectangular region in Fig. [6] (b) Magni- 
fication of rectangular region in Fig.[7fa). The fractal structure of the 
chaotic repeller can be clearly seen. 

where q is a real-valued parameter. The scaling behaviour of 
the partition function 

r(e, q) ~ eb-WM (18) 

defines the generalized or Renyi dimensions d(q). The dimen- 
sion d(q) can be determined numerically for given q ^ 1 by 
computing the partition function T as in Eq. (T7\ for increas- 
ing box sizes e M = 2~ M with M = 10,9,8,7,6,5. The 
weight of each resulting box is then the sum of the weights of 
the corresponding smaller boxes. Finally, the dimension d is 
extracted by fitting the function ( fT8l . The dimension d(Q) is 
equal to the Minkowski dimension, also called capacity or box 
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counting dimension. In the limit q — > 1, the Renyi dimension 
becomes the information dimension 

= lim 7— w k i Inwki ■ (19) 
fez 

Similar to the general case of with (77^ 1, can be de- 
termined by fitting d(l) lnejvf to ^ fc; wu In as function 
of the box size em = 2~ M . The dimension d(2) is called cor- 
relation dimension since it is related to correlation functions 
as, e.g., in Ref. 15711 . 

In the case of a fractal set with uniform measure, i.e., a 
point belongs to the set (/ = 1) or not (/ = 0), all dimensions 
d(q) have the same value. In the more general case of fractal 
sets with nonuniform measure (/ e R + ) the dimensions d(q) 
differ. In such a case one speaks about a multifractal. The 
chaotic repeller is strictly speaking a set with uniform mea- 
sure. For numerical reasons, however, we have seen that it is 
more convenient to consider a real-valued I(s,p). That this 
approximation is justified can be seen by the fact that the nu- 
merically computed fractal dimensions d(Q), d(l), and d(2) 
are all approximately 1.68. The predicted exponent according 
to Eq. (f2]l and d = 1.68 is then around 1.84. Unexpectedly, 
this is not close to the Weyl exponent a € [1.96,2.02] ob- 
tained from the counting of resonances. Hence, it seems that 
the conjectured fractal Weyl law fails for the dielectric mi- 
crostadium. 

In order to see the origin of this failure, let us consider 
the escape rate 7 of the chaotic repeller. Points in the vicin- 
ity of the repeller escape exponentially in time with rate 7. 
To determine this rate, we remark that the numerically com- 
puted repeller in Fig. [6] is an approximation to the real re- 
peller, in that numerically computed phase space points are 
never exactly located on the repeller. We can exploit this 
fact to determine 7 by evolving the numerically computed 
repeller in time. Its spatial structure does not change un- 
der time evolution, only the total intensity decays exponen- 
tially. The corresponding decay rate is nothing else than the 
escape rate 7. We find 7 0.056c/i?. This translates into 
Im(f2) = —jR/(2c) k —0.028 which is significantly larger 
than the mean escape rate of the quasi-bound modes in Fig. [4] 
We conclude that the chaotic repeller as defined above also 
fails to describe the mean lifetimes of the quasi-bound modes. 

In the following we show that both problems stem from 
the fact that the conventional chaotic repeller ignores the par- 
tial escape at dielectric interfaces according to Fresnel's laws. 
Rays entering the leaky region are not completely transmitted 
to the exterior of the cavity. For TM polarization the reflection 
coefficient is 

R™=(^-^) 2 (20) 
Vsm (x + Xt)J 

with the angle of incidence \ an d the angle of transmission \t 
related by Snell's law n sin \ — sin \ t . We can include Fres- 
nel's laws in the concept of the chaotic repeller by allowing 
for a real-valued intensity / assigned to each ray. A similar 
approach has been used to describe far-field emission pattern 
based on unstable manifolds j48l l55ll . Initially, we choose 



1=1 for a given ray, and whenever the ray enters the leaky 
region at a phase space point (s,p) we multiply the intensity I 
by the reflection coefficient Rtm(s,p)- Finally, the ray's in- 
tensity is < I < 1. The chaotic repeller defined in this way 
is a set with nonuniform measure I(s,p); it is a multifractal. 

Using the extended concept of the chaotic repeller we com- 
pute the escape rate to be 7 w 0.034c/i?. This gives Im(f2) « 
—0.017 which is in much better agreement with the distribu- 
tion of imaginary parts of the quasi-bound modes depicted in 
Fig- El Figure[8]shows the chaotic repeller including Fresnel's 
laws. The fractal dimensions turn out to be d(Q) ks 1.986, 
d(l) w 1.913, and d(2) » 1.877. The predicted exponent for 
the Weyl law is therefore 1.99, 1.96, and 1.94, correspond- 
ingly. The values resulting from the box-counting dimension 
and the information dimension are in the interval of the Weyl 
exponent a G [1.96, 2.02] calculated from the resonances. For 
the high-index microstadium we therefore have confirmed the 
conjectured fractal Weyl law using the extended concept of 
the chaotic repeller. 
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FIG. 8: Chaotic repeller of the microstadium with n = 3.3 including 
Fresnel's laws; cf. Fig. [6] Intensity is higher for darker regions, and 
vanishes in the white regions. 



V. LOW-INDEX MICROSTADIUM 

Now we proceed to the microstadium with low refractive 
index n = 1.5. Such a cavity has been fabricated recently by 
using a PMMA polymer matrix j47l l58ll . Here, we can com- 
pute the resonances in the interval Re(fi) £ [0, 75] as can be 
observed in Fig. [9] For example, RefO) = 75 corresponds 
to R = 7.2 /im for A = 600 nm 1471 15811 . The comparison 
with Fig. [3] reveals that the resonances of the low-index sta- 
dium are much deeper in the complex plane, which reflects 
the stronger degree of openness due to the larger leaky region 
I sin xd — V n - This can also be seen from the statistics of 
imaginary parts depicted in Fig.[l0] Here, we observe a clear 
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Re(£2) FIG. 11: Number of modes TV with cutoff C = 0.3 versus frequency 

Q in a double-logarithmic plot (solid line) for the microstadium with 
FIG. 9: Resonances in the complex-frequency plane for the microsta- n = 1.5. The slope of the linear fit is a w 1.75 (dashed line; offset 
dium with n = 1.5. In total 2893 resonances have been found. for clarity). 
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FIG. 10: Probability density of imaginary part of resonant frequen- 
cies £1 for the low-index microstadium. 



Figure QT| shows a fit to the Weyl law ([TJ with C = 0.3 
yielding a « 1.75. Varying the cutoff parameter C in the 
interval [0.25, 0.4] we find a E [1.68, 1.88]. For the strongly 
open microstadium with low index of refraction we see a clear 
deviation from the conventional Weyl law for closed systems 
as one would expect. 

The chaotic repeller including Fresnel's laws is shown in 
Fig. Q~2] The escape rate for phase space points near the 
chaotic repeller is determined to be 7 w 0.26c/ R. This trans- 
lates into Im(£7) w —0.13 in good agreement with the distri- 
bution of the quasi-bound modes in Fig. [10] Note that Fres- 
nel's laws do not change much here since the reflectivity in the 
leaky region is in general small due to the low refractive in- 
dex. The numerically computed dimensions are d(0) » 1.512 
and g?(1) w g?(2) w 1.593. Hence, the predicted exponent is 
1.76 and 1.78 which is fully consistent with the fractal Weyl 
law with exponent a G [1.68, 1.88]. 
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FIG. 12: Chaotic repeller of the low-index microstadium including 
Fresnel's laws. 



VI. CONCLUSION 

We have demonstrated that the harmonic inversion tech- 
nique allows to compute very efficiently a huge amount of 
complex resonances in dielectric microcavities. This is of 
high practical value for the statistical analysis of resonances in 
the context of quantum chaos and complex scattering. The ap- 
proach has been illustrated for the chaotic microstadium made 
of GaAs and polymers. We demonstrate that the fractal Weyl 
law can be applied to the class of dielectric cavities by extend- 
ing the conventional concept of the chaotic repeller by includ- 
ing partial escape of light according to Fresnel's laws. We ex- 
pect that our finding is of high relevance for the understanding 
of ray-wave correspondence in chaotic microcavities. More- 
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over, it can have implications on multi-mode lasing in such research group 760 is acknowledged, 
cavities Il59ll . 
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